Numerical studies of the extended two-chain model of friction 
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We investigate numerically a simple microscopic model to describe wearless dry friction between 
atomically flat contact interfaces without thermal fluctuations (T = OK). The analysis of the 
incommensurate ground state shows a breaking of analyticity when the amplitude e of the Lennard- 
Joncs interaction between the both harmonic chains increased beyond a critical value e c . By the 
introduction of a suitable order parameter and using a finite size scaling we could show that the 
breaking of analyticity is a second order phase transition. This transition is often called Aubry 
pH | transition. By applying of an uniform external force F ex we have determined the critical force of 

. static friction F c above a sliding motion occurs. 
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^ , I. INTRODUCTION 
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J~ There are several ways for studying friction phenomena. Experimental observations are very powerful to charac- 
i | terize the frictional behaviour and to formulate empirical friction laws, e.g., Coulomb- Amonton's laws |Q, from a 
phenomenological point of view. Macroscopically, everybody can usually observe that during two solids are slided 
against each other, wear take place. But wear is not necessary for the occurrence of friction. In the last few years the 
technological development makes it possible to study wearless dry friction between atomically flat contact interfaces 
experimentally (nanotribology) |l]-|3| . These nano-scale experiments on systems with less complicity stimulated again 
the theoretical efforts to investigate friction phenomena on the base of microscopic models H-jl^l- The dissipative 
nature of friction is a typical non-equilibrium problem. From the theoretical point of view, wearless dry friction be- 
tween two atomic flat solids could be explained by using a simple driven mechanical many-body system with a lattice 
structure and take into account interaction forces. Thus the question arises, how the structure of such a microscopic 
model of friction has to look. 
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t^J- . One of the simplest structure on the microscopic level is given by the discrete version of the well-known Frenkel- 
Kontorova (FK) model |l3| , which is in general a model of an one-dimensional harmonic chain of particles in a periodic 
potential. The FK model is mainly used for describing an adsorbate atom monolayer on an atomically flat clean 
substrate surface. Using this picture, the driven FK model, where a constant external force is applied to each particle 
of the chain, can also be used as a simple model structure for a wearless friction situation, which has been studied by 
several works ■ Note, several physical phenomena such as the dislocations motion in crystals |L4j , commensurate- 
incommensurate transitions (l^Jl^l and charge-density waves (CDW's) |lT|]l8fl could also be investigated by reducing 
on the FK model structure. Definitely, the simple FK model is a good starting point in order to construct more complex 
and realistic friction models on the atomic level. Thus, a natural extension of the FK model is the replacement of the 
rigid substrat (hard body) by a deformable substrat monolayer pinned on a bulk (soft body). Such a more complex 
O ■ model of friction was firstly introduced by Matsukawa and Fukuyama [0J^]. They considered an one-dimensional 
model consisting of two deformable interacting harmonic chains, where each particle of the lower chain is harmonically 
• *h , pinned on a rigid bulk. The static and kinetic properties of this two-chain model of friction was analysed in great 
details by Kawaguchi and Matsukawa |) 10 1. So far, the FK model and the two-chain model are simple models of 
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wearless friction between an atomic monolayer and a hard or soft body, respectively. For describing a more realistic 
friction situation, which should represent wearless friction between atomically flat bodies, it is necessary to replace 
the adsorbate monolayer by a substrat. For this reason, Weiss and Elmer (l^Jl^l proposed an interesting approach, 
the so-called Frenkel-Kontorova-Tomlinson (FKT) model, which is an one-dimensional lattice model for a soft upper 
body sliding on a hard lower body. In contrast to the two-chain model, with the FKT model they extended the FK 
model structure in another direction by harmonically coupling of each particle of the adsorbat monolayer on an upper 
macroscopic sliding mass, whose position relative to the lower surface is characterized by the coordinate xb- Note, 
all of these models are investigated as completely one-dimensional microscopic friction models, where the particles of 
the one-dimensional chain can only move parallel to the direction of the external force. Now, it is obviously clear, 
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how you have to extend the previous model structures in order to write down a more complex microscopic model 
system, which is necessary to describe friction phenomena as simple but also realistic as possible. For that purpose 
we introduce an extended two-chain model consisting of two interacting harmonic atomic chains embedded in a two- 
dimensional space, where each particle of the upper or lower chain is harmonically pinned on an upper or a lower 
rigid bulk, respectively. This model describes two soft bodies with an one-dimensional contact interface, which can 
carry out a genuine two-dimensional sliding motion against each other. In our studies we consider the extended two 
chain model as a simple relaxation model by taking energy dissipation via a phenomenological damping term in the 
mechanical equation of motion into account. Hence, all model, where the energy dissipation is modeled by hand 
through an explicit damping term, cannot explain the origin of friction in complex physical systems refering to the 
sliding state. Well, that is also not the aim of our investigations, because with our approach we want to establish the 
typical macroscopic friction situation (tension experiment) on a microscopic level. Nevertheless, the introduced model 
structure can be used to discuss the mechanism of the occurrence of static friction, which is one characteristic feature 
of dry friction, based on an energy concept. The hope is that the properties of the ground state - defined as the 
minimal potential energy state of the mechanical model system in absence of an external force - describe completely 
the static friction behaviour as function of the interaction between the two deformable chains. In general the ground 
state structure can exhibit two situations as a driven force increased from zero. On the one hand, a zero static friction 
arises exactly then if the maximum static friction force vanished, i.e., the model system is frictionless. On the other 
hand, a finite static frition corresponds always with the existence of a finite maximum static friction force, which is 
the maximum potential energy state of the system above which a sliding state appears. In other words, if we evolve 
the particle configuration from their ground state configuration, then the maximum static friction force is defined by 
the largest depinning force, which is equivalent to the smallest driven force above which no stable stationary state 
exists. 

For the analysis of static properties of the extended two chain model we follow the philosophy of Aubry by using 
strangely enough an incommensurate ground state structure, which has the only purpose to introduce a special 
kind of mathematical functions called as hull- functions. On the basis of the specific hull-function concept for the 
visualisation of the ground state structure of the discrete FK-model in the incommensurate case Aubry yJHJ has 
been shown that a new type of phase transition (called as breaking of analyticity) occurs when the amplitude of the 
periodic potential increased beyond a critical value. That means, below a critical amplitude the hull-functions have a 
continuous (analytic) structure whereas above a critical amplitude the hull- functions are discrete (non analytic). In 
our numerical investigations of the hull- functions of the extended two chain model we have also found the transition 
by breaking of analyticity (Aubry transition). In order to distinguish, whether the Aubry transition is really about a 
second order phase transition, we have specifically studied this type of phase transition like Coppersmith and Fisher 
(but which have been investigate the behaviour near the maximum static friction force jl7],[l8|]) from the viewpoint of 
critical phenomena. For that purpose we have introduced a suitable order parameter, which shows the typical critical 
behaviour for second order phase transitions and critical exponents can be calculated via scaling arguments. 

Here, it should mentioned, that from the physical intuition a finite maximum static friction as function of the 
interaction strength e has to exist in all previous introduced models equal whether one choose a commensurate or an 
incommensurate ground state structure. The main differences are the following. In the commensurate case the static 
friction is always larger than zero for any e > 0. By way of contrast, in the incommensurate case whether the static 
friction is equal or unequal zero depends strongly on e below or above a critical value e c due to the preculiar breaking 
of analyticity. Hence the commensurate and the incommensurate ground state structure reveals two types of physical 
system classes characterizing by different critical behaviour in the sliding state near the threshold force F c . The first 
class is given by model systems with a commensurate ground state structurce. Effectively, commensurate models are 
characterized by finitely many internal degrees of freedom, which lead to trivial critical exponents in the sliding state 
for the velocity-force characteristic near the threshold F c . An incommensurate ground state structure is equivalent 
to a second class of models with infintely many internal degrees of freedom. Well, both model classes emphasize 
the fact that the nature of static and kinetic friction is deeply rooted with internal degrees of freedom generated by 
a complex atomic interaction behaviour, which is the key for unterstanding dissipation in general. Commensurate 
model have one main disadvantage, because an uniquely visualisation concept for the ground state structure like 
the hull-function for the incommensurate case is missing. Moreover, commensurate models exist always in a pinned 
ground state structure, that means, models imply static friction for any value of the amplitude of the interaction 
strength larger then zero intrinsicly. Therefore, for understanding the physical fundamental mechanism of the arising 
of static friction it seems to be most interesting to study incommensurate models, because the ground state shows 
two different incommensurate phases related to amazing properties. That means, an unpinned incommensurate phase 
corresponds to zero static friction whereas the pinned incommensurate phase is equivalent to a finite static friction. 

In the present paper we introduce the extended two chain model and the numerical calculation methods. The 
detailed investigations of the static properties, especially the ground state structure of the extended two chain model 
are discussed by defining an order parameter and proposing scaling relations. 
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II. THE MODEL SYSTEM 



Following Matsukawa and Fukuyama @,^|, we introduce the extended two-chain model, which is a two-dimensional 
microscopic lattice model system describing wearless dry friction between two atomically flat solids. Therefore, we 
consider two deformable atomic chains, so called an upper and a lower chain. The atoms in each chain are classical 
point particles and are able to realize a two-dimensional movement. We assumed a harmonic particle interaction in 
each chain. The interaction between the two chains is of the Lennard-Jones type. Moreover, both chains are coupled 
harmonically on an upper or a lower substrat respectively. Like Ref. j7j the effects of energy dissipation in the upper 
and lower chain are phenomenologically approached to be proportional to the difference between the velocity of the 
i-th particle and that of the center of gravity. The external force F_ ex — (F ex , 0) acts on the upper chain in x direction 
without y component. In order to simplify the situation, we consider overdamped motion. The equations of motion 
for the particles in the upper and the lower chain are given for each direction as follows (see also FIG. I) 

N b 

m ala ( x? - (if) ) = fc* [ xU + < +1 - 2<] - A£ 6 [< - zc a ] + ]T Ff nt « - x b ) + F™ (I) 

3=1 

m ala ( y? - (#> )=kl[ yU + V?+i - 2tf] - Kit, [Vt ~ (h + L Q )]+J2 Kt (tf ~ Vj) ( 2 ) 

3=1 

N a 

m blb ( x\ - (x\) )=k* b [ xU + x b +1 - 2x\] - k b / ub [x\ - zc b ] + £ *Zt {A ~ 4) ( 3 ) 

3 = 1 

m b l b ( y\ - (yl) ) = kf[ y\_ x + y\ +1 - 2y\) - k b s l b [y\ - b b ] + £ t (y b - y") , (4) 

3 = 1 

where rf = (xf,yf) and r b — (x b ,y b ) are the immediately position vectors of the i-th particle in the upper and 
lower chain. Furthermore, N a (N b ), m a (m b ) is the particle number and the particle mass, 7 a (75) denotes the friction 
coefficient and c a (c b ) is the mean atomic spacing in ^-direction of the upper (lower) chain. It should be mentioned, that 
the particle mass and the friction coefficient are for all particles in each chain the same. For a compact representation 
one can introduce the chain index I — (a, b) relating to the upper and lower chain. The velocity of the center of 

gravity is given by the average C^i') = 5jEil r i-^/i where X- stands for x\ and y\ respectively, fcf and k\ are 
the strength of the harmonic interaction force in the chains for the x and y direction. The strenght of the coupling 
of each chain on an upper and a lower substrat is represented by k^ b and fc^, which determine the rigiditiy of the 
chains. A further system parameter is the substrat-substrat distance L = b a + b b + Lq, where b a (b b ) is the mean 
atomic spacings in y-direction between the upper (lower) chain and upper (lower) substrat. Lq denotes the distance 
between the two harmonical chains and is a function of their interaction strength. The interaction force between the 
upper and lower chain, projected on the x and y axes, can be written as 

E^K-<)=E^)^^ (5) 

3=1 3=1 13 



Nj ( t 



^Ktbd-vj) =E^3-)^-^ (6) 



3=1 



with the chain indices / = (a, 6), J = (b, a) and the distance = \rf — r b \ = J (x b — x") 2 + (y b — y^) 2 relating to 
the z-th atom of the upper (a) and the j-th atom of the lower (b) chain. F(rij) can be expressed by the potential 
V(rij) through F(rij) — — d V qJ 13 ^ ■ As interaction potential we choose specifically the Lennard-Jones potential 

representing the character of the friction due to the intermolecular forces in the best way. Here e is the strength 
of the interaction and a is a characteristic length parameter, which parameterized the minimum of the potential 
ymin _ V(rij)\ r .. =ra via the distance from equilibrium r = a2 1 ^ e . 
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III. NUMERICAL METHOD 



For solving the two-dimensional extended two-chain lattice model we choose periodic boundary conditions in both 
chains 

x \ = Xi+Nj - L I ( 8 ) 

Vi = Vi+Ni ( 9 ) 

with the chain index I = a,b and the chain lengths Li = Njci. This implies that the system size of the upper and 
lower chain are the same L = N a c a = N^cj,. Hence, the ratio or the misfit of the two lattice constants could be 
expressed by the particle numbers 

c a N b 

a= — = —. 10 

Cb N a 

This means, that the ratio a of the two lattice constants is always a rational number in the numerical realisation. 
Using the concept of the hull-function, which was introduced by Aubry and coworkers |l9| , |20|| for describing the ground 
state of the underlying model system, we need an irrational ratio a between c a and q,. If we have a rational a, i.e., 
in the most commensurable case a = 1, we cannot apply the concept of the hull- function. Because the whole ground 
state configuration of each chain, consisted of Nj particle positions, would be represented by one value of the argument 
of the hull function and hence the hull-function plot lose his mathematical power of proposition. For this reason, it is 
necessary to investigate an incommensurate ground state structure relating to the ratio a of the two lattice constants. 
Unfortunately, in the numerical calculation we could only approximate an irrational number through optimal rational 
ratios of two natural numbers using a finite continued-fraction expansion up to a certain order |lq , [l9| 

1 , . 

a = a Q H j (11) 

a\ H 1 



a 2 



1 

a n -x H 



Following the literature, we choose the simplest continued-fraction expansion with ao — 1 and a\ = a 2 = • • • = a n = 1, 
so that a could be expressed by the ratios of the Fibonacci numbers F n+ i = F n + F n -i, where n = 1, 2, 3, • • • as well 
as Fq = 1 and F% = 1. Hence, 

F n+2 _ 3 5 8 13 21 34 55 89 144 233 377 610 987 1597 
Q ~ F n+1 _ 2'3'5'Y'1^'2T'34'^'W'W4'2^'3T7'610' "987 ' " ' ^ ' 



and the irrational limit a = linin^oo = i (\/5 + l) of the approximation series (12) is the inverse of the so-called 
golden mean. From the point of number theory Jl9| , this is the most incommensurable case of the underlying model 
system, which could be emulated in a systematical way. In our studies we set the lattice constant of the lower chain 
equal one (c& = 1). Accordingly, the upper mean lattice spacing c a is equal a and we calculated with the system sizes 

N b 89 233 377 610 1597 



N a 55 ' 144 ' 233 ' 377 ' 987 



(13) 



It should be mentioned, that not all ratios oft the approximants are "good" irrational number approximations in the 
same manner because for certain particle number ratios there exist some artifacts in the hull function. 
For the numerical calculation of the equations of motion we used a fourth-order Runge-Kutta algorithm with an 
adaptive stepsize control. In order to discuss the essential features of the present model we restrict the wide range 
of model parameters by setting m; = 1, 7j = 1, kf = krf = 1, fc^ b = = 1 and a = 1. As initial conditions 
we assumed at t = a particle configuration in the x-y-plane, where the iV a and Nb chain particles are arranged 
at regular lattice sites periodically. Besides, the distance Lq is set to an arbitrary initial value, i.e. Lq = 1. For 
calculating the stable stationary states we used the same criterion like Kawaguchi and Matsukawa (To) . That means, 
we used a velocitiy condition 



\ ' {N a + N b ) 



< 10~ 5 (14) 
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for stopping the Runge-Kutta calculation. Hence, the system could be considered to be in a static state. The next 
step consists in the check up of the total sum of the coupling forces in the (a)-chain and (b)-chain in x-direction 

N a N b 

F2 b = Y.KZA^-ic a ] , ^ = 51 -3<*] (15) 

i=i j=i 

and y-direction 

N b 

(h + L )} , F b s Z b = Y,k^ b [y]-b b ] (16) 



p Qy _ ,a y r a _ 

r sub — 2-^1 K sub 111% 
1=1 



respectively. In absence of a driven force (F ex = 0) it is sufficient to consider only the value of the total sum of the 
coupling forces ( |l6| ) in y-direction, because for t — 0, after the Runge-Kutta steps and in the final state the ^-direction 
forces ( fig) are equal to zero. The regular periodical arrangement of the particles for the initial time t = leads to 
Fgub = F S ub = 0- After the relaxation in the final stationary state relating to the velocity condition (|l4|) the total 
sum of coupling forces (|l6| ) must always be zero, which is in general not true for an arbitrary chosen initial value 
Lq. Accordingly, during the numerical calculation every final stable state in respect to Lq must be proofed and if 
necessary Lq must be changed like Lo = Lq ± 5Lq. Then the chains relax again and the procedure would be repeated 
until the system reached a final stationary state, which fulfil the condititon F^ b = F s * b < 5 • 1CT 5 (—0). In this 
way we could find numerically the stationary configuration of the extended two-chain model and for F ex = one can 
assumed that this is approximately the ground state, defined as having an energy which cannot be lowered by moving 
of some system particles. The properties of the ground state could be analysed by changing the interaction strength 
e, which corresponds to set the correct Lq — Lo(e) numerically. 

In the presence of an external force F ex > acting in ^-direction one must slightly change the numerical calculation 
procedure. Instead of F^ b = F^ b = the x-direction forces have for t = the values F^ b = F ex and F^ b = 



respectively. After each Runge-Kutta step we have to guarantee that the control parameters (15|) have again these 
values. Hence, we are able to evolve the model system from the ground state into the last possible static solution by 
continuously increasing of the applied external force. In this sense we could numerically define the threshold, where no 
longer a static solution of the underlying system exist. This threshold F ex = F c is the critical force of static friction, 
which is also called maximum static frictional force. Below F ex < F c the chains are in a steady state, such that the 
afterwards space-time averaged velocity v = (if), t — 0. Above F ex > F c the upper chain, which is after an initial 
transient decay in a steady sliding state, slides with a constant velocity v. Note, in an easy way one can investigate 
the present model in a gravitational field via the constraints (p"6|). 



IV. NUMERICAL RESULTS 



A. Ground state 



First, we investigate the stationary states of the extended two-chain model in the absence of an external force, i.e., 
pex _ q -pjjg stationary states are given by the solution of equation system (|l]-[|) with vanishing right-hand side 

= m a7a (r? - (if) ) = m blb (rf - <r^) ) (17) 

with consideration of the conditions (|l5|) and ([r]), where F^ b = F^ b = F^* b — F^ b = 0. Hence, the ground state of 
the extended two chain model is the stationary particle configuration having the minimal potential energy for F ex = 0. 
The knowledge of the ground state determine a physical system, such as the static friction, completely. 



1. The hull-function 



Similarly to the FK model jl9| or FKT model |Tl[] , the lattice structure of the particle positions in the incommensurate 
ground state can be uniformly described by a hull-function. For examining the lattice structure of the present two- 
dimensional extended two-chain model, it is necessary to introduce four hull-function due to two degrees of freedom 
of the particles in each chain. The two hull-functions h a and h c of the upper chain are defined as 
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xl = ic a + (f> + h a (ic a + (p) (18) 

yl = b b + L + (f ) ' + h c {i c a + (p 1 ) (19) 

and the two hull-function h b and h c of the lower chain are given by 

x\ = i c b + + h b (i c b + i/j) (20) 

V* = b b + ij' + h d (ic b + ip') , (21) 

where <\>, cf>' , ip and ip' are constant phases. For our investigations, it is conviencent to set cf) = </)' = tp = ip' = 0. The 
hull-functions are periodic and even, which can usually be expressed by 

K(z + Cb) = h a (z) = —h a (—z), h c (z + Cb) = h c (z) = —h c (-z) (22) 

h b (z + c a ) = h b (z) = -h b {-z), h d (z + c a ) = h d (z) = -h d (-z) . (23) 



The argument z of h a and h c is defined in the range from to Cb, whereas of h c and h d it is given from to c a - Note, 
the visualisation concept of the hull-function makes it necessary that the argument of the just defined hull-functions 
h a , h c or hb, h d must be the same, respectively. Because, the displacements in x- resp. y-direction of all particle 
in each chain can be assigned to different values of the argument z due to the irrational ratio a of the mean lattice 
spacing c a and c^- 

First, we have numerically calculated the ground state of the extended two-chain model for a fixed ratio a — Nb/N a 
of the particle numbers given by N a = 233 and Nb — 377. Several investigations of the FK model [^0|, FKT model 
p"T[ such as the two-chain model jl(| showed, that the ground state strongly depends on the interaction strength e. 
Therefore, we started the computation for different values of e in the wide range from 0.05 to 0.4. For small values 
of e up to a certain value e c we get the result that the periodic hull- functions h a , hb, h c and h d are analytic, i.e, 
smooth and continuous. This behaviour is illustrated in Fig. 2, where the interaction strength has the value e = 0.22. 
Furthermore, the stationary state for e = 0.22 is characterized by the distance Lq = 0.986273 between the upper and 
lower chain. Figure 3 shows another numerical calculation of the hull-functions for the interaction strength e = 0.33 
with the corresponding distance L$ = 0.959584. We can see that the hull-functions in Fig. 3 are discrete, which 
means that they are no longer analytic. This discrete structure is related to jumps, which occur at certain points of 
the argument z of the hull-functions. Now, the comparison of Fig. 2 with Fig. 3 shows the phenomenon of breaking of 
analyticity, which is due to two significant aspects of the hull- functions h a , hb, h c and h d for the interaction strength 
e < e c and e > e c , which defined a certain transition value e c above the hull-functions are no longer analytic. For a 
detailed analysis of the threshold e c we introduced a suitable order parameter in the following part. 

Well, Aubry has shown (20) that the breaking of analyticity of the incommensurate ground state of the FK model 
is characterized by the existence of a largest central gap at the half of the period of the hull-function. Whereas the 
situation in the two-chain model of Matsukawa and Fukuyama jl(J is quite different for several elastic parameters and 
three simple cases are observed: both hull- functions h a and hb have a largest gap at the half of the period (FK limes), 
whether h a or hb has only a largest central gap and finally both hull-functions do not have a largest central gap but 
two symmetrical gaps with the same size refering to the half of the period. In Fig. 3 we can also observe that no 
largest central gap for all hull-functions of the extended two-chain model exists. Furthermore, the hull-functions show 
also two symmetrical gaps of the same size refering to the half of the period of the hull-function in the whole analysed 
parameter range. Up to now we have only investigated the ground state structure for a fixed system size (N a = 233 
and Nb = 377) and a certain range of interaction strength e. Now, we vary the system size of the underlying model 
which corresponds the different approximation orders of the irrational number through a rational ratio of the particle 
numbers. Here, we calculated especially the ground state structure for the ratios a = 89/55, 233/144, 610/377, and 
1597/987 in the same range of the interaction strength e like above. In particular, we proofed whether there occurs an 
anormalous behaviour in the structure of the hull-functions for certain values of e due to the using of a rational ratio 
instead of an irrational number in the numerical computation. Indeed, Fig. 4 shows such a anormalous behaviour of 
the hull-functions for e — 0.33 and for the particle numbers N a — 377 and 610. In the comparison with Fig. 3 we can 
recognize that the even particle number Nb = 610 generates a wrong particle position at the half of the period x = 0.5 
of hb in Fig. 4 due to the periodical boundary condition. Consequently, in term of a suitable rational particle ratio 
Fig. 3 shows the right picture of the hull-functions for e = 0.33, whereas Fig. 4 shows some artifacts in the hull-function 
structure for the same interaction strength e = 0.33 as a result of bad choice N a = 377 and N b = 610 of the system 
size for the interaction strength e = 0.33. Beyond it, there appear much more artifacts in the hull-function structure 
different ranges of e als consequence of a bad choice for the particle number ratios. For N a = 144 and N b = 233 the 
occurence of such artifacts destroyed the order parameter analysis of the transition area. Kawaguchi and Matsukawa 
interpreted in [ |To[ these artifacts, which also occur in the two-chain model for certain values of the interchain strength 
and ^V a = 144, N b = 233, as information of the appearance of a new gap structure at the elastic parametes Kb- This 



6 



interpretation is wrong, because, if they calculated the model system with the same interaction strength for different 
particle numbers these artifacts will be vanished. The appearence of artifacts in the hull-function for the ground state 
defines the range of validity of the hull-function concept for several e in the numerical calculation, which makes it 
necessary to proof the hull-function structure for the same e for several system size. Relating to stationary state which 
is given by Eq. (|l7|) , Fig. 5 shows the distance Lq between the upper and the lower chain as a function of the strength 
e of the interaction force, which was calculated for the particle numbers N a = 233 and Nb = 377. The plot in Fig. 5(a) 
shows the computation of Lq for e in range from 0.05 to 0.45 and Fig. 5(b) shows Lq in the transition area. Note, the 
distance Lq is remarkably insensitiv to the system size. For the system sizes a = 377/233, 610/377, and 1597/987 
we found a deviation for Lq — Lo(a) smaller than 10~ 5 . In the occurence of local artifacts in the hull- functions for 
certain particle numbers and certain interaction strength the deviation for Lq(o) lies in the order 10~ 4 . Note, that in 
the ground state above the transition point e c also many local rearrangements take place which can be identified as 
first order phase transitions. 



2. The order parameter 

Now, we analyse the transition of breaking the analyticity near the threshold e c in terms of second order phase 
transition. For this it is convenient to introduce an order parameter, which is an indicator for the occurence of a 
phase transition at a critical point. The order parameter distinguishs both phases which are involved in the phase 
transition, because it is defined to be zero in one phase. Refering to the transition of breaking of analyticity the order 
parameter is defined in the discrete phase, which measures the discrete jumps in the hull- functions. Therefore the 
order parameter is in the continuous phase for e < e c equal to zero and in the discrete phase above e c unlike zero. 
Here, we use the definition 



OP m = 



d m f(z) 



dz 1 



dz 



5(0) 



(24) 



Here, 5(0) behaves like 5(0) ~ Az 1 |az^o • If the hull-function is smooth in the interval [0, 1], the integral becomes 
finite and the order parameter vanished due to the prefactor 6(0) _1 . On the other hand, if f(z) has a discrete 
step at 00) w e obtain d z f(z) ~ 6(z — z a ) and therefore J Q [d z f(z)] 2 dz ~ J 6(z — zq fdz ~ 6(0)J Q 6(z - z ) ~ 6(0), 
i.e. the order parameter becomes finite. Due to the finite system size, which devides the region of integration in 
equidistant subintervals Az — l/N a (Az = 1/Nb) for the upper (lower) chain, the integral in Eq.(p4|) reduced to a 
sum over the particle numbers. Therefore, we have intrinsic discontinuities in the numerical computation, because the 
hull-functions are given by a discrete set of points at Az at all, which define the systematical error of the calculation. 
In order to smeared out the systematical intrinsic discontinuities, we choose for our numerical investigation OP4 as 
order parameter. Hence, we get the following expression 



OP = OP A = 



d i f(z) 



dz 4 



dz 



5(0) 



(Azf 



-Na-2 

£ 

k=2 



( I fk+2 - 4 A+i + 6 f k - 4/ fc _i + h-2 I ) : 



(25) 



It should mentioned, that OP shows necessary for all hull-functions h a , h c , h c and ha the same qualitativ behaviour. 
Therefore, in order to show all basic features of the phase transition (transition of breaking of analyticity) it is enough 
to calculate OP (via Eq.(E51)) only for the hull- function h a = h a (z). 



3. The scaling behaviour 

Next we have calculated the order parameter given by Eq.(p5|) for different system sizes a = 89/55, 233/377, 610/377 
and 1597/987. Figure 6(a) shows the order parameter as function of the interaction strength e in the transition region. 
We clearly see that the order parameter OP is changed from zero below a certain critical threshold e c to a value unlike 
zero above e c for all particle number ratios. This is a certain hint that the transition of breaking of analyticity is a 



second order phase transition. Using a finite-size scaling [E3 24 for the data in Fig. 6(a) given by the scaling function 



OP(e-e c ,L) = L*OP[L7(e-e c )) (26) 
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we find that the order parameter OP exhibits a scaling behaviour in the transition region shown in Fig. 6(b). That 
means, all curves of Fig. 6(a) could mastered on one resulting curve (Fig. 6(b)). The system size L is determined 
by the particle numbers 55, 233, 377 and 987. In the transition region, which can be characterized by the reduced 
interaction strength t = (e — e c )/e c , we find that the data of Fig. 6(b) are consistent with the following power law for 

e > e c 

OP(e-e c )~^( £ -e c f (27) 

with the critical exponents f3 = 0.239 ± 0.05, S = 0.687 ± 0.05 defining the ratio fj, = /3/S w 0.347 ± 0.05. Due to 
the finite-size scaling we obtain for the threshold interaction strength of the transition of breaking of analyticity the 
critical value e c = 0.25575 ± 0.001. 

B. The static friction 

In the following we want to calculate the maximum static friction force F c of the extended two-chain model refering to 
a pinning structure which is given in the incommensurate case for e > e c . It is clear that in the presence of an applied 
forces F ex > ground states are no longer defined. From the physical intuition, we hope that one can adiabatically 
evolve the underlying model system from the ground state by increasing F ex from zero, because the threshold force 
F c is the largest depinning force F ex above no a longer stationary state exist, i.e., Eq.(|l7|) is no longer valid. Like 
Ref. [^0| the static friction is given by 

N a N b 

^fl = -EE( F ** « " 4) ) t = N *< PBX )t (28) 

which can be obtained from the equation of motion (^j)-(^) by summing over i and averaging with respect to t. It 
should remarked, that Eq.(|2^) is also valid for the calculation of the kinetic friction in a sliding state, i.e., for F ex > F c . 
For the determination of the threshold force F c , we choose especially the interaction force e = 0.33 and the system 
size N a = 233 and Nb = 377. Note, the numerical calculations of F c for the extended two-chain model are very hard 
and the computation time is very large. In a first calculation, we find that the maximum static friction force has 
approximately the value F c 0.4 but the numerical deviation is very large. Therefore, it takes much more efforts in 
order to calculate the accurate value of the threshold force. 



V. CONCLUSIONS 

We have investigated the ground state for a two-dimensional extended two-chain model with an incommensurate 
lattice structure. We have found that the structure of the hull-functions depends strongly on the strength e of the 
interaction between the upper and the lower body. It was numerically observed the breaking of analyticity in the 
ground state for the critical value e c = 0.25575. Furthermore, it was shown that this type of phase transition could 
be uniform described by an order parameter. This order parameter representation allows a systematic analysis of the 
critical behaviour for e > e c due to the calculation of the critical exponents. Consequently, the peculiar breaking of 
analyticity shows all features of a second order phase transition. 

A second interesting result of the ground state structure occurs in associtation with the numerical computation 
for different interaction strengths as function of the rational system sizes. It has been shown that artifacts in the 
hull- functions structure appear for certain rational particle numbers of a fixed interaction strength e. This fact pointed 
out that the correct ground state structure depends strongly on the choose of sufficient rational ratio of the particle 
numbers as a irrational number approximation for a certain range of the interaction force. 

Concerning the static properties of the extended two-chain model in the presence of an applied forces F ex > we 
have calculated the maximum static friction force F c . 

Finally, it should remarked, that the extended two-chain model is a simple model structure , which established the 
macroscopic dry friction situation on the microscopic level between atomically flat surfaces by using a simple driven 
mechanical many-body system. It is convenient to define the model as simple as possible in order to dicuss some 
basic features like the ground state structure or the maximum static friction behaviour. Therefore we took in our 
investigation a simple relaxation dynamic into account and neglected thermal fluctuation (T = K) . Surely, such 
kind of microscopic models could be always extended by additional thermal fluctuations representing by stochastical 
forces which leads to creep in the presence of applied forces. Note, the vality of the underlying model system is given 
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by the assumption of simple harmonic particle interaction. It is clear that this assumption in many real systems where 
large distortions occur. 

Note, there are one possible and very interesting extention of the extended two-chain model by the introduction 
of quenched randomess which corresponds surface impurities. It is an open question, whether randomess destroy the 
frictionless state in the incommensurate case for e < e c in the extended two-chain model like in the FK- model 
Well, it is an interesting question, whether randomess, which exists in any real system, leads in general to a finite 
maximum static friction force for any interaction strength e. Refering to the compare whith experiments, another 
natural extention of the present model is to proceed to a three-dimensional system, which is able to describe wearless 
dry friction between two atomically flat plans pinning in each case on a substrat and embeding in a three-dimensional 
space. 

In later works we want to investigate the threshold behaviour near F c especially for F ex > F c . It should be proofed 
whether the velocity-force characteristics shows a critical behaviour or not. Well, for a complete understanding of 
the behaviour near the threshold F c a renormalization group calculation is needed, but surely very complicated for a 
deterministic mechanical system. 

Note in the present paper, we have considered a simple overdamped motion. Hence, it is also interesting to discuss 
the influence of inertia terms (underdamped dynamic) on the behaviour near the threshold F c and on the sliding 
state. 
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FIGURES 



FIG. 1. A schematic picture of the two-dimensional extended two-chain model. The external force F ex is given by the 
applied force F Fex which acts on each particle of the upper chain times the particle number N a of the upper chain 
where N a , (Nb) is the particle number of upper (lower) chain. 

FIG. 2. The hull- functions of the ground state plotted for the particle numbers N a = 233 and N b — 377 with the 
interaction strength e = 0.22. For an uniform representation all hull-functions are renormalized by the mean lat- 
tice spacing c b and c a , respectively. Therefore, the substitution for the hull- function h a and h c of the (a)-chain is 
z —* z' = z I 'cb, whereas for the hull-functions hb and hd of the (b)-chain it is z — > z' = z/c a . 

FIG. 3. Hull-functions of the ground state for e = 0.33 and the particle numbers N a = 233 and Nb = 377. 

FIG. 4. Artifacts in the hull-functions for a bad choice of the particle numbers, e.g., given by iV a = 377 and Nb = 610 
and the interaction strength e = 0.33. 

FIG. 5. The interaction distance Lq in y-direction between the upper and lower chain as function of the interaction 
strength e. The distance L is plotted in different ranges (a) from e = 0.05, • • • , 0.4 and (b) from e = 0.24, • • • , 0.31 
(transition range). 

FIG. 6. In the computation it is convenient to rescale the order parameter by OP — > OP - (Az) 6 . The order parameter 
OP = OP(e — e c ) describes the critical behaviour near the threshold e c . The graph (a) shows the behaviour of the 
order parameter for different particel numbers and (b) is the finite-size scaling plot of the mastered order parameter, 
where OP(e)L~z is plotted versus (e — e c )L>i . 
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